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Abstract 

In this study, we investigate hybrid utilization of angle-of-arrival (AOA) and received signal strength 
(RSS) information in visible light communication (VLC) systems for 3D localization. We show that AOA- 
based localization method allows the receiver to locate itself via a least squares estimator by exploiting 
the directionality of light-emitting diodes (LEDs). We then prove that when the RSS information is 
taken into account, the positioning accuracy of AOA-based localization can be improved further using a 
weighted least squares solution. On the other hand, when the radiation patterns of LEDs are explicitly 
considered in the estimation, RSS-based localization yields highly accurate results. In order to deal with 
the system of non-linear equations for RSS-based localization, we develop an analytical learning rule 
based on the Newton-Raphson method. The non-convex structure is addressed by initializing the learning 
rule based on 1) location estimates, and 2) a newly developed method, which we refer as random report 
and cluster algorithm. As a benchmark, we also derive analytical expression of the Cramer-Rao lower 
bound (CRLB) for RSS-based localization, which captures any deployment scenario positioning in 3D 
geometry. Finally, we demonstrate the effectiveness of the proposed solutions for a wide range of LED 
characteristics and orientations through extensive computer simulations. 
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I. Introduction 

The reliable positioning in the environments where the Global Positioning System (GPS) signals cannot 
easily penetrate has become a growing research area due to its numerous applications I!]. Among many 
other technologies, radio frequency identification (RFID) |j2l-|3], ultra wideband (UWB) positioning 
Il5ll . cellular or wireless local area network (WLAN) based positioning IQ, Q, and fingerprint based 
localization techniques iH are the most prominent ones for GPS-free localization. Beside those tech¬ 
nologies, visible light communication (VLC) l!9l|- lfTTII . which has been recently emerging as a promising 
technology to provide high data rates for short range communications, can also deliver high accuracy for 
3D localization ifT^ . IT^ . 

VLC systems have numerous advantages in comparison to other legacy systems. First, the performance 
of VLC systems is not affected by the radio frequency (RF) interference. Moreover, multipath fading can 
be averaged out as the area of photodetectors are very large compared to wavelength of the visible light 
llT4l . Second, light emitting diodes (LEDs) offer narrow beamwidth, which enables more accurate angle 
of arrival information at the receiver side ifTSl . lIThll . Third, VLC systems can be utilized in scenarios 
where RF radiation is potentially hazardous or even forbidden such as planes or hospitals. Lastly, VLC 
systems maintain maintains high data rate communications while simultaneously providing illumination. 

Due to its aforementioned advantages, employing VLC systems for the localization is recently getting 
significant attention in the literature. In particular, two approaches have become prominent for VLC-based 
localization: angle of arrival (AOA) based methods |[T5l . |[T6l and received signal strength (RSS) based 
methods ifTTl . llT8]l . While AOA-based approach takes the direction information of LED transmitters, 
RSS-based localization considers the captured power from EED transmitters. Eor example, in ITSl . AOA- 
based localization is investigated for indoor scenarios and a connectivity-based solution is proposed under 
2D settings. The proposed solution exploits the narrow field of view (EOV)s of EEDs and relies on that 
fact that the receiver is likely to be located at the intersection of the directions of the connected EEDs 
transmitters. Yet, the proposed method is valid only under 2D settings where two straight lines always 
intersect between each other unless they are not parallel to each other. In ifThl . AOA information of 
EED transmitters are exploited in an ad-hoc networking environment. By starting from a limited number 
GPS-enabled devices, each node discovers its own location and broadcasts this information in order to 
help the other nodes to find their own locations. In ifTTl . an RSS-based indoor positioning algorithm 
based on preinstalled EED transmitters is proposed to estimate receiver locations by analytically solving 
the equations that characterize the Eambertian pattern lfT9l . However, the analytical derivations are not 
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generic, and cover only a limited number of scenarios, where the directionality of the LEDs (i.e., their 
modes) and the receiver’s direction are fixed to some predetermined values. In lITSl . Cramer-Rao lower 
bound (CRLB) for RSS-based localization is derived for a specific scenario of LED fransmitters and 
receiver positions, where the authors emphasize the complexity of obtaining the CRLB in a 3D geometry. 

In this study, we evaluate AOA-based and RSS-based localization methods by using multi-element 
visible light access points (VAPs) which consist of several LED transmitters ifTOl . If20l . ll^ . The main 
contributions of this paper are as follows: 

• Hybrid Localization: We develop hybrid localization methods that utilize both AOA and RSS 
information in the estimation process to improve the positioning accuracy in VLC systems. Lor 
AOA-based localization, we show that the VLC receiver can estimate its own location via a least 
squares (LS) estimator and the positioning accuracy can be improved further by incorporating with 
the RSS information. In addition, we demonstrate that the outcome of AOA-based location estimation 
can effectively be useful to deal with the non-convex objective function in RSS-based localization. 

• Theoretical Lramework: We introduce a comprehensive theoretical framework that allows us to 
investigate the various hybrid localization methods for VLC and the impact of the orientations 
of LED transmitters and the physical characteristics of LEDs on the estimator performance. As 
opposed to earlier work in the literature, the proposed methods and theoretical investigations are 
generic, and applicable to any 3D topology. Based on the introduced framework, we discuss the 
optimal LED weighting problem for AOA-based localization and the non-convex objective function 
for RSS-based localization. As a benchmark, we derive the analytical expression of CRLB for 
RSS-based localization, which generalizes the derivations provided in ifTSl . By using the derived 
CRLB, we demonstrate the trade-offs between light distribution and localization accuracy distribution 
and discuss the impact of LED characteristics, such as directivity, density, and orientation, on the 
maximum achievable positioning accuracy. 

• Learning Algorithms: We introduce two algorithms, i.e., analytical learning algorithm based on 
Newton-Raphson method and random report and cluster (RRC) algorithm, to combat with the 
system of non-linear equations and the non-convex structure of the objective function for RSS- 
based localization, respectively. In order to increase the likelihood that the Newton-Raphson method 
converges to the global optimum, we employ RRC algorithm and show its effectiveness numerically. 

The rest of the paper is organized as follows. The system model which captures any VAP deployment 
scenario in 3D geometry is provided in Section ini AOA-based localization and hybrid utilization of AOA 
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Fig. 1. System model with multiple VAPs and a single VLC receiver. 


and RSS information for VLC systems are discussed in Section HIT] Subsequently, RSS-based localization 
is investigated in Section |IVl Theoretical expression of CRLB for RSS-based localization, the analytical 
learning rule, and the RRC algorithm are provided in this section. The numerical results evaluating the 
performance of the proposed methods for different configuration of VAPs are given in Section |Vl Finally, 
some concluding remarks are indicated in Section |Vll 

Notations: Ik is the K x K identity matrix. The transpose operation is denoted by (•)^. The Moore- 
Penrose pseudoinverse operation is denoted by The operation of 11-112 is the 2-norm of its argument. 
The vectorization and the expectation operators are denoted by vec (•) and E [•], respectively. The kernel 
of a matrix is represented by ker (•). The field of real numbers is shown as M. M (0, C) is fhe Gaussian 
distribution with zero mean and covariance matrix of C. 


II. System Model 

We consider K VAPs communicating with a VLC receiver as illustrated in Fig. [T] Without loss of 
generality, the location of the receiver and its orientation are denoted by fr = [x,y,z]^ G and 

Dr = G in Cartesian coordinate system, respectively. We consider multi-element 

VAPs which consist of M LED transmitters. The location of mth LED transmitter of feth VAP and its 
orientation are denoted by = [xrak^Vrak^ZjnkV ^ and S 

respectively. By assuming transmit power of 1 W for each EED transmitter, the signal power of the mth 
EED of kth VAP at the receiver is given by IT9l 

Pmk COs(0niA:) 


27r 


p2 

X^mk 


\0fovJ V ^/2 


( 1 ) 


August 25, 2015 


DRAFT 
















5 


where (pmk is the angle between the LED transmitter orientation vector and the incidence vector, Omk is 
the angle between receiver orientation vector and the incidence vector, Rmk is the distance between the 
LED transmitter and the receiver, Ar is the area of photo detector (PD) in m^, 0 fov is the EOV of PD, 
n is the mode number which specificities the directionality of LED, and If (•) is the rectangle function 
defined as 


n(x) = < 


for |x| < 1 
for |x| > 1 


( 2 ) 


While n (^mfc/^FOv) in 0 implies fhat a VLC receiver can defect the light only when O^k is less than 
0FOV^ n ((/>mfc/(7r/2)) ensures that the location of VLC receiver is in the EOV of LED transmitter. Let 
Vmfc = tr — Tmk = CmkV G be the incidence vector between the receiver and the mth 

LED transmitter of A:th VAP. Then, the parameters of ([1]) can be expressed as 



Rmk — ||'^mA:||2 ; 


( 3 ) 

COS(^772/c) 


'^mk^mk 

( 4 ) 


Rmk 

11V mk 112 

COs(^6mf^^ 

_ n^ir mk - tr) _ 
Rmk 

II2 

( 5 ) 


Therefore ([T]) can be rewritten as 

{jl + 1)Ar / A TT 4’mk\ f f \ 

=-(W j " Ij " 

where 


r . (vIfcn^fc)”v^fcnR 

= {^mknUf^ + brnklT-Uj^ + Cmk^Uf^) 

T 

^ r”9 79 n ' n+3 • to 

i^mk + bik+(^mk) = 

It is worth noting that ® includes only the line-of-sight (LOS) component, i.e., first order Lambertian 
pattern, since the LOS component of signal is dominant compared to other multipath components for 
VLC systems ifTTl . lfT9l . Hence, we discuss the theoretical analyses given in the following sections based 
on LOS component. In this study, we assume that all LEDs in VAPs are identical. In addition to that, the 
locations of the transmitters and directions of each LED are assumed to be known at the receiver, which 
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Fig. 2. The distance between VLC receiver’s location and the line extended by the direction of an LED transmitter. 

can be shared through periodic broadcast messages. The unique identity of each LED transmitter and 
its corresponding VAP are assumed to be decodable at the receiver, which can be achieved by assigning 
different codes/headers from each LED transmitter. Accordingly, it is assumed that the receiver is able 
to measure the RSS associated with each LED transmitter. 

III. AOA-Based Localization 

In AOA-based localization, a VLC receiver selects one of the LED transmitters for each VAP based on 
the RSS information and locates itself by exploiting the orientations of LED transmitters. Geometrically, 
AOA-based localization relies on finding a point in 3D space such that it minimizes the sum of distances 
(or squared distances) between this point and all the other lines extended from the directions of LED 
transmitters selected by the receiver ifTSl . ifThl . This approach can allow the receiver to locate itself with 
only two anchor nodes when the receiver is close to the intersection of the corresponding LED directions. 

In order to analyze the AOA-based localization, let be the line associated with the mth LED 
transmitter of /cth VAP as illustrated in Pig. |2l where its origin and direction are captured by the vectors 
Tjnk and rimfc, respectively. In addition, let be the projection matrix which projects every 

vector X € onto the null space of n^k, i-e., ker(nmfc)> which can be calculated by A^k = 

I 3 — UrnkT^mk- Hcnce, any vector in the column space of Amk is orthogonal to the direction of L^k- In 
particular, the projection of any point on L^ak yields the same point, i.e., the intersection between L^k 
and the subspace spanned by the columns of A^^. The intersection point h^k can be simply calculated 
as Amk^mk sincc Tmk is a point on L^k- Hence, the distance vector between an arbitrary point 0 € 
and is given by = ^mk - A^^Q, which yields 

Ayyil^Q -|- . ( 8 ) 
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By stacking the set of equations related to LED transmitters, we obtain 

b = A0 + d, (9) 

where 


bsii 




dsii 

bsfcfc 

, A = 

Asfcfc 

, d = 

dsfefc 

PskK_ 


J^skK_ 


_^skK_ 


and Sk is the index of selected LED transmitter of kth VAP. Therefore the objective function which 
minimizes the sum of squared error is solved via an LS estimator given by 

fR = Atb. (11) 


A. Weighting LED Transmitters 


Treating the AOA information obtained from each of the LED transmitters identically, as in (fTTI) . may 

cause significant positioning error since the distance between exact position of the receiver and the line 

extended by the direction of LED transmitter may be large in some cases. In other words, less reliable 

AOA information may bias the location estimation and degrade localization accuracy. In order to address 

this issue, we consider an objective function which minimizes the weighted sum of squared distance as 

K 

fR = argmin ^ I3mk\\dmk\\l 
k=l 

m=Sk 

K 

= argmm ^ /3mfc||A^fc0 - b^fclll , (12) 

k=l 

m=Sk 

where pmk is the weighing factor for the distance between 0 and Lmk- The purpose of weighing factors 
is to take the variability of into account in the optimization. Under the assumptions of E [d^fc] = 0 
and E [d^^d^fc] = 0, according to Gauss-Markov theorem ll22l . the optimal weighting factors that yield 
minimum variance estimator for (fT 2 l) is obtained as 


l^mk Cq X 


1 


(13) 


]E[||d^fcP] • 

where cq is an arbitrary positive real number. In order to calculate (fTSl) . the statistical characteristics of 
dmfc should be known a priori. To this end, we exploit RSS information associated with LED transmitters 
in this study. 


August 25, 2015 


DRAFT 









Without loss of generality, consider an LED transmitter located at the origin and its orientation is 
set to no = [0,0,1]^. Assuming that the LED transmitter and the VLC receiver face each other, i.e., 
Hr = — no, 0FOV = 7’'/2, and < vr/2, we rewrite ® as 


Po = 


(n + 1 )Ar 


^n+l 


(14) 


271" + 2 /^ + 2 

We then obtain the identity which characterizes the distance between the LED direction and VLC 
receiver’s location as 


dolll = C'"+3 


Z n+3 — z 


(15) 


where ||do ||2 = and C = ■ H is worth noting that (fTSl ) also identifies the locations 

where VLC receiver could observe the same RSS and forms an RSS contour. Lor example, when Pq 
is fixed to 0.1 mW and Ar = 1 cm^, the potential locations of a VLC receiver are shown in Lig. |3 
for a given LED mode. Under the assumption that the receiver is equally likely to be on the locations 
that satisfy (fTSl) . one should obtain E [UdoHl] in order to find opfimal weighfs given in ([T3l) . However, 
fhe calculafion of E [UdoHl] is nof frivial due to nonlinear characteristics of (fTSl) . On the other hand, 
considering the narrow beamwidth of the LED with high modes, RSS contours are spread around on 
z-axis as shown in Lig. [3] Therefore, it becomes more likely that the VLC receiver location far away 
from the LED direction. Hence, in this study, we approximate the exact distribution of ||do ||2 with a 
uniform probability density function (PDL) given by 

1 


Pdidolb) = < 


27Tdu 

0 , 


||do||2 — da 
otherwise 


where 


dmax = max (lido II 2 ) = V B 


-C-B- 


-c 


(16) 


(17) 


and B = (2n + 6)/(2n + 2). Therefore the variance of ||do ||2 is obtained as approximately. By 


induction, the optimum weights in (fT3]) are derived as 

2-KPmk 

Pmk 


Co X 


{n+l)A-RZ 


— = Pmk 


(18) 


B-— - B- 

where (a) follows from the appropriate selection of cq. As a result of (IT^ . LED transmitters are weighted 
based on their RSS information Pmk- 

The problem given in (O corresponds to an unconstrained quadratic optimization problem, which can 
be solved via LS method. The closed-form solution of (IT^ can then be obtained as 

tr — A^v^b^y , (19) 
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z [m] 


Fig. 3. RSS contours for a given LED mode. All the locations on the RSS contours have identical RSS values. 


K K 

where ^ ^ Pmk-^mk und ^ ^ Pmk^mk- 

k=l fc=l 

m=Sk m=Sk 

IV. RSS-Based Localization 

Under the assumption of availability of physical characteristics of LEDs, the receiver can locate its own 
location based on RSS information associated with LED transmitters. Eet s G be the observation 

vector which is given by 

s = p(0) + n , (20) 

where n G ^ ^ ^2j km) is an additive Gaussian noise vector, p(0) = vec (P(0)) G 

and P(0) G R^x^ is a matrix which contains the exact RSS information with the entry on mth row and 
/cth column being Pmk- The log-likelihood function for the location of VEC receiver is then expressed 
as 

/:(0) = logp(s;0) , (21) 

where 

^ 2^®^^ (^-^(S-P(0 ))^(s-P(0))) (22) 


dPmk 

dCLrnk 

dPmk 

dbmk 

dP-mk 

dCynk 


(n+l)AR 

ni 

1 ^mk 

27r 

^FOV 

(n+l)AR 

n( 

^mk 

2 it 

^ 0FOV 

(n+l)AR 

n( 

^ Omk 

Stt 

> ^FOV 



(29) 
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and 0 G is the parameter vector which corresponds to the location of the VLC receiver, i.e., fr. 
Maximum likelihood (ML) estimate of fr is therefore formulated as 

fr = argmax£ (0) oc argmax (—(s — p(0))"'"(s — p(0))) . (23) 

As a result, (|2^ can he expressed as a nonlinear least squares (NLS) problem given hy 

fr = argmin (||s - p( 0 )|| 2 ) . (24) 

0 

In the following subsections, first, we introduce a learning rule in order to solve (l24l) . We then investigate 
the non-convex structure of (l24b and address the initialization issues of the learning rule via the RRC 
algorithm. Finally, we provide CRLB of (l24b . 


A. Learning Rule for RSS-Based Localization 

In order to solve the system of nonlinear equations in (l24b . we follow multivariate Newton-Raphson 
method, which yields 

gi+i ^ 0 i _ (s - p( 0 *)) , (25) 


where ry G (0,1] is the step size and H is the Jacobian matrix of p(0) with respect to 0. As 0i, 
03 corresponds to x, y, and z in (l20l ). respectively, H can be explicitly given by 



dPoo 

dPoo 

dPoo 

dx 

dy 

dz 

dPio 

dPio 

dPio 

dx 

dy 

dz 


dPfAK 

dPfAK 

OPmk 

dx 

dy 

dz 


02 ^ and 


(26) 


As it can be seen in (l26l) . each row of H indicates how the RSS associated with an LED transmitter 
change when the receiver moves in one of the axes, i.e., x, y, and z. Considering the chain rule of 
derivatives, the row associated with the mth LED transmitter of kth VAP is calculated as 



Since v^fc = fr —the matrix in (|27]) . i.e., Jacobian of respect to 0, becomes an identity matrix. 
Therefore, (|27]) can be directly calculated by evaluating the derivative of Pmk with respect to at 
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Fig. 4. Convexity of in. In this setup, there are three regions where the modeling error is less than le-11 and lead to three 
different local optimum points which indicate non-convex structure of J24l l. While the first initial point allows l l25t to find the 
global optimum, i.e., receiver’s location, the second initial point converges to one of the local optima. 


Vmfc = tr — r^fc- The derivative of P^k with respect to is analytically given in (l29l) hy employing 
an auxiliary function g (x; a, b, k, I, m, re) defined as 


g (x; a, b, k, I, rre, re) = 


df (x; a, b, k, I, rre, re) 
dx 


{ax + k)^^~^\c3X^ + C2X^ + CiX + Cq) 

(x^ + rre)^ 

where / (x; a, b, k, I, rre, re) is a function based on the structure of Pmk and keeping one of the elements 
of incidence vector v^fc as a variable and the others as constants in ([7]). It is explicitly defined as 

(ax + k)^'^~^\bx + 1) 


f (x; a, b, k, I, rre, re) = 


(28) 


(x^ + rre) 2 

where cq = kbm + lamn, ci = bmna + bma — lk{n + 2>), C 2 = —{kb{n + 2) + 3la), and C 3 = —2ab. Affer 
plii§§iii§ ^mk ^mk into (|29|), i.c., a^nk X ^mk^ ^mk U Umk^ and 21 ^mk^ ^261) and 


(l27l) are obfained fheorefically. Therefore, the learning rule given in (1251 ) can be calculated analytically, 
which eliminates the numerical calculation of H and yields less complex structures. 


B. Non-Convex Structure of RSS-Based Localization 

Although Lambertian pattern offers a convex set by itself, (|2^ is not a convex function in general. 
This is due to the fact that the set of feasible solutions associated with each LED transmitter may become 
closer to each other in multiple locations in 3D geometry. For instance, consider a3mx3mx3m 
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where the VAPs are located at the corners as illustrated in Fig. |4l In this setup, we assume that each VAP 
has single LED transmitter where n = 30 and the angles between LED directions and side walls are set to 
45 degrees. The VLC receiver parameters are set as 0 fov = 85 degrees, Ar = 1 cm^, fr = [1,1,0.75]"'", 
and riR = [0,0,1]^. Having individual RSS information from each LED transmitter based on ([T]), we then 
find the locations where VLC receiver could observe the same RSS information for each LED transmitter 
in 3D geometry, which corresponds to RSS contours. 

As it can been seen in Lig. |4l the VLC receiver is located at the intersection of RSS contours. However, 
there are other locations in the room where the RSS contours get closer to each other. When we evaluate 
(l24l) and find fhe locations where fhe modeling error, i.e., ||s — p(0)||2, is less fhan le-11, we observe 
fhree differenl regions which indicate at least three different local optimum points. Therefore, different 
initial points for the learning rule given in (|25] ) may converge to different location estimates. Lor example, 
in Lig. m while the first initial point yields the true location of the VLC receiver, the second one converges 
to one of the local optima, i.e., the point where the RSS contours get closer to each other. 

C. Initialization of RSS-Based Localization 

In order to increase the likelihood that (1251 ) converges to the global optimum, one may use the solution 
of AOA-based localization method as an initial point. The AOA-based solution is useful, such that it 
provides a closed-form solution that is close to the exact location of the VLC receiver, and hence does 
not require exhaustive search methods. However, this method may still lead to wrong results in some 
cases. In order to identify better initial points which increase the likelihood of finding the global optimum, 
i.e., location of the receiver, we propose a heuristic algorithm which is based on random search If23l and 
unsupervised clustering methods |[24l . which we call as the RRC algorithm. 

RRC algorithm includes two steps: 1) random reporting and 2) clustering. In the first step, the algorithm 
generates S random location samples in 3D space. The rewards achieved by the location samples are 
then calculated based on the modeling error in (|2T]) . By evaluating the rewards, only G < S points that 
give the most promising results are chosen and the rest of the sample locations are eliminated. If the 
solution space is sampled in a sufficiently dense way, the favorable location samples should be located 
in the space as clusters since the objective function is non-convex and local regions appear in 3D space 
as clusters as illustrated in Lig. |4| 

In the second step, we aim at finding fhe cenfroids of fhe clusters. The reason for finding fhe cenfroids 
is fhat fhey are more likely to be around the true location of the VLC receiver, which make them good 
candidates for the initial points. To this end, we exploit unsupervised learning methods which are well 
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investigated in the machine learning literature. Due to its simplicity, in this study, we employ /C-means 
algorithm which is an effective method for finding the centroids of the clusters in a high-dimensional 
space |[2^ - lf26l . Essentially, K-means algorithm converges to the centroids of the clusters with iterations. 
It applies two steps per each iteration. In the first step of jth iteration of if-means algorithm, the points are 
partitioned into C clusters, , R\= 2 i ■ ■ ■ ’ based on their distances to the centroids obtained 

in {j — l)th iteration. In the second step, the fth centroid, i.e., c*, is updated to better fit for the location 
samples as 


„0'+i) 


1 



’ 


(29) 


where pi is the £th point obtained from random reporting step and | • | is the cardinality of its argument. 
Once we obtain the centroids, we utilize them as initial points for the update rule in (1251) . We then select 
the best point based on (|2^ . 


D. Cramer-Rao Lower Bound for RSS-Based Localization 

Without loss of generality, when x ~ A( (/r(0), C(0)), the element on ith row and jth colunm of 
Fisher Information Matrix (FIM) J(0) can be calculated as |[22l 


[J(0)].i = 


\dp{e)] 

T 

C(0)“^ 

'Opiey 

[ 09, \ 


[ 09, \ 


1 




dOi 


de, 


(30) 


Since the noise term in (l20l) is assumed to be white Gaussian noise and does not depend on the location 
of receiver, the second term is zero and first term yields [J(0)]ij = ^ ■ Based on (l20l) . 

/r(0) corresponds to mean power captured power from all transmitters, i.e., /r(0) = p(0), which leads 
to the following expression: 


1 

r9p(0)i 

T 

■5p(0)' 


L 09, \ 


L oej \ 




By using (iMl) . complete FIM can be then calculated as 

J(0) , 


(31) 


(32) 


where H is the Jacobian matrix of p(0) with respect to 0, which is already derived analytically by using 
(l27l) and (l29l) . Once J(0) is calculated, CRFB for the receiver location in 3D geometry is obtained as 


RMSE{x, y,z) > y/tr(J-H0)) 


(33) 


August 25, 2015 


DRAFT 



























14 



Fig. 5. Simulation setup. 

where RMSE(x, y, z) is the root mean square error (RMSE) of the estimator given in (|2^ . It is worth 
noting that the analytic expression of H allows arbitrary scenarios in 3D coordinates, which generalizes 
the theoretical results provided in llT8l where the angle between detector plane and floor plane is fixed 
to 0 degrees. 


V. Numerical Results 

In this section, we evaluate the performance of AOA and RSS based methods and impact of physical 
characteristics of LED transmitters on the estimators through computer simulations. Eor simulation 
tractability, we consider an empty room where its depth, width, and height are set to 5, 4, and 3 meters, 
respectively. We deploy 4 VAPs which are located at the comers of the room as shown in Eig. [5] The 
angle between the direction of VAP and ceiling is denoted by ^ceiling- Each VAP has 4 LED transmitters 
where the angle between the direction of VAP and the direction of LED transmitters is assumed to be 
identical and parameterized as 0poiar- For the VLC receiver, 0 fov and Ar are set to 85 degrees and 
1 cm^, respectively. The direction of the VLC receiver is assumed to be riR = [0,0,1]^. Einally, based 
on the noise model given in iflSl and the references therein, we set to le—13 A?. 

A. Performance of AOA and RSS Based Location Estimators 

In order to evaluate the estimators discussed in Section |III] and Section |IV] numerically, we consider 
two different paths, i.e.. Path 1 and Path 2, as shown in Pig. |5] Path 1 considers a VLC receiver where 
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(a) Path 1, n = 10. (b) Path 1, n = 30. 




(c) Path 2, n = 10. (d) Path 2, n = 30. 


Fig. 6. Positioning accuracy on Path 1 and Path 2 (0poiar = 20°, ^ceiling = 30°). 


its localization on x-axis and y-axis are fixed to 2 m. On the other hand, Path 2 follows a horizontal path 
where the exact position of the VLC receiver on y-axis and z-axis are set to 1 m and 1.5 m, respectively. 
We then sweep the VLC receiver’s position on z-axis for Path 1 and x-axis for Path 2 and calculate 
RMSE of each estimator after 50 realizations. We perform the simulations for n = 10 and n = 30 in ([T]) 

when ^ceiling — 30 and 0polar — 20 . 

As it can he seen from Fig. 0 LED mode, i.e., n, does not affect the performance of pure AOA method 
based on (fTTl) for both Path 1 and Path 2. This result is expected since the LED transmitters are assumed 
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to be identical. Therefore VLC receiver selects the same LED transmitters as anchors regardless of n, 
which yields the same estimation results for pure AOA method. However, LED directivity, characterized 
by re, significantly affects the performance of AOA-based estimation, when RSS information is taken 
into account. Eor weighted AOA method based on ([T^ . increasing re provides better estimation results 
on both Path 1 and Path 2. This is due to the fact that RSS information becomes a dominant factor when 
VLC receiver is close to the line pointed by the orientation vectors of LED transmitters. In other words, 
nodes with higher re have more reliable AOA estimates for a given RSS value, which is captured by the 
proposed weighted LS estimator. 

Lor RSS-based localization, we consider two different initializations for the Newton-Raphson method 
(r/ = 0.2) given in (1251) . We first run the simulations when weighted AOA results are considered as the 
initialization points. In this case, the search algorithm may happen to converge to local optima, which 
causes significant positioning errors. Lor example, the estimator performance degrades drastically for 
Path 2 when the VLC receiver is close to the side walls of the room, as given in Lig. | ^b)[ On the other 
hand, when RRC algorithm (5 = 500, G = 100, and C = 4) is applied for finding betfer initial points, 
RSS based localization is more likely to find the global optimum and attains the CRLB in most of the 
positions. 

When the position of VLC receiver is close to the ceiling, all methods that use RSS fail as shown 
in Lig. [ ^a)| and Lig. | ^b)[ This is due to the large angle between receiver orientation vector and the 
incidence vector which reduces the effective area of PD. In particular, after certain height, LOV of VLC 
receiver does not allow the PD to capture any signal, which explains why CRLB goes to infinity. 

B. Number Of Clusters for RRC Algorithm 

In this subsection, we evaluate the impact of number cluster of RRC algorithm, based on the configu¬ 
ration given in Lig. |5] As a pessimistic scenario, we consider a VLC receiver located at the origin of the 
coordinate system. We then provide the probability of convergence to the location of VLC receiver for 
C = {0,1, 2,3,4} when re = 10 and re = 30. Lor each scenario, we also include the result of weighted 
AOA to the initial points obtained via RRC algorithm. In order to have a better understanding on the 
convergence, we evaluate the impact of number clusters in the noiseless case. We assume that a successful 
converge is considered when Hfr — frUI is less than le—2. 

As shown in Lig. |7j increasing number of clusters yields better probability of converge. However, 
different LED modes require different number of clusters. Lor example, while RRC with 2 clusters is 
sufficient to attain the global optimum with very high probability when re = 10, RRC with 4 cluster. 
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Fig. 7. Probability of convergence for RRC algorithm. 

i.e., 4 extra initial points beside weighted AOA result, yields high probability of convergence to global 
optimum when n = 30. 

C. Configuration of VLC Access Points 

In this subsection, we investigate the impact of deployment of VAP on positioning accuracy. We 
consider the simulation setup illustrated in Fig. [5] and choose ^ceiling as a design parameter. We sweep 
^ceiling between 20° and 70° and calculate CRLBs for the positions of the VLC receiver based on 3D 
grid in the whole room. We then calculate the probabilities where CRLB is less than or equal to certain 
values, i.e., CRLBq = {0.01,0.04,0.07,0.125,0.25} in meters when n = 10 and n = 30. In the analysis, 
we fix 0polar to 20°. 

As it can be seen from both Fig. l g{a)| and Fig. l g{b)| the optimum ^ceiling that gives the highest 
probability of accuracy in the room varies depending on CRLBq and LED mode n. For example, the 
highest probability of accuracy is obtained at ^ceiling = 60° when CRLBq = 0.01 m as given in Fig. I 31a)[ 
On the other hand, the highest probability is achieved at ^ceiling = 30° for the same LED mode when 
CRLBq = 0.25 m. Similarly, when n is set to 30, the choice of ^ceiling = 37° yields the highest probability 
of accuracy for the same CRLBq as given in Fig. | g{b)[ 

It is also possible to infer how LED directivity affects the positioning accuracy overall in the room 
from Eig. [8] As it can be seen from Eig. [ g{a)[ %90 of the room is covered when CRLBq = 0.25 m and 
n = 10. However, increasing LED mode n to 30 reduces the probability of accuracy to %85 of the room 
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(a) n = 10. 



(b) n = 30. 

Fig. 8. Probability of accuracy for a given 6^ceiiing (^poiar = 20°). 


as given in Fig. | g|^b)[ On the other hand, higher LED mode increases the prohahility of highly accurate 
results. For example, when CRLBq is fixed to 0.04 m, %45 of the room is covered for n = 30 while it 
reduces to %27 of the room for n = 10. 

D. Orientations of LED Transmitters 

In this subsection, we discuss the impact of the orientations of LED transmitters on positioning accuracy. 
Considering the setup illustrated in Fig. |5j we sweep 0poiar between 5° and 40°. Similar to the analysis 
given in Section IV-Cl we obtain the probabilities that CRLB is less than or equal to certain values when 
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(a) n = 10. 



(b) n = 30. 

Fig. 9. Probability of accuracy for a given 6'poiar (^ceiling = 30°). 


n = 10 and n = 30. For this analysis, we fix ^ceiling to 30°. 

As shown in Fig. |9l the optimum 6*poiar depends on the LED mode and selected CRLBq. For instance, 
when n = 10 and CRLBq = 0.01 m, the optimum 0poiar is 17.5° as given in Fig. [ 2Ia)| Flowever, as 
shown in Fig. [ 2th)[ the optimum 0poiar is obtained as 10° when n = 30. Similarly, when CRLBq is set 
to 0.25 m, 27° and 19° yield to the highest prohahility of accuracies in the room when n = 10 and 
n = 30, respectively. Nevertheless, in general, the optimum 0poiar increases for higher CRLBq. This 
is due to the fact that larger 0poiar yields more physical separation between LED transmitters for each 
VAP and the intensity distribution becomes more homogeneous in the room. Therefore while narrow 
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0poiar yields more accurate results at specific points of the room, larger 0poiar provides less accurate but 
better probability of accuracy for the setup given in Fig. |5] On the other hand, increasing 0poiar above 
the optimal point deteriorates the positioning accuracy as the orientations of LED transmitters start to 
be parallel to the walls of the room. In these cases, the center of the room does not receive sufficient 
amount of energy, which causes more positioning error. 

VI. Conclusion 

In this study, we discuss AOA-based localization and RSS-based localization methods by considering 
hybrid utilization of AOA and RSS information in the location estimation. We show that AOA-based 
localization can be solved with an LS estimator. Yet its estimation results may be highly inaccurate, 
depending on the location of VLC receiver. On the other hand, when RSS information is utilized to 
weight LED transmitters in the optimization, the positioning accuracy increases significantly. Eor the 
scenario investigated in this study, i.e. a room with the dimensions of 3 x 5 x 4 meters and 4 VAPs, the 
localization error is less than 1 meter. On the other hand, RSS-based localization method, which exploits 
the Lambertian patterns of LEDs, offers high positioning accuracy at the expense of a system of nonlinear 
equations and a non-convex objective function. In order to solve the system of nonlinear equations, we 
derive an analytical learning rule based on the Newton-Raphson method. Since the learning rule is 
analytical, it eliminates the calculation of Jacobian matrix numerically. Lor the hybrid utilization, the 
learning rule is initialized with the result of AOA-based method to increase the likelihood of converging 
to the global optimum, which the positioning accuracy improves up to 10 cm. In addition to employing 
the results of AOA-based method, we also develop a heuristic search method, RRC algorithm, to identify 
extra initial points which could lead to find the global optimum. 

In this investigation, we also discuss the impact of the orientations of LED transmitters and the physical 
characteristics of LEDs on the localization performance probabilistically. Lor this purpose, we utilize 
the CRLB that is derived for an arbitrary configuration in 3D geometry in this study. According to our 
analyses, when the illumination is homogeneous in 3D geometry, the positioning error becomes relatively 
high but mostly in attainable levels. On the other hand, when the illumination is not homogeneous due to 
the orientations of LEDs or using highly directive LEDs, the positioning accuracy is improved significantly 
at the locations where the energy is high; this also degrades the positioning accuracy at the same time 
in the locations where the energy, i.e., illumination, is low. 
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